Weak antilocalization in high-mobility two-dimensional systems 
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Theory of weak antilocalization is developed for high-mobility two-dimensional systems. Spin- 
orbit interaction of Rashba and Dresselhaus types is taken into account. Anomalous magnetoresis- 
tance is calculated in the whole range of classically weak magnetic fields and for arbitrary strength of 
spin-orbit splitting. Obtained expressions are valid for both ballistic and diffusive regimes of weak 
localization. Proposed theory includes both backscattering and nonbackscattering contributions 
to the conductivity, ft is shown that magnetic field dependence of conductivity in high-mobility 
structures is not described by earlier theories. 
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I. INTRODUCTION 

Anomalous magnetoresistance caused by weak local- 
ization is a powerful tool for extracting kinetic and band- 
structure parameters of three-dimensional (3D) and 2D 
systems.^ Theoretical expression for magnetoconductiv- 
ity valid in the whole range of classically weak magnetic 
fields taking into account all interference processes has 
been first derived in Refi. In the absence of spin-orbit 
interaction the sign of the magnetoresistance is negative. 

However the anomalous magnetoresistance is an alter- 
nating function in 2D semiconductor systems. In partic- 
ular, in low fields it is positive and cannot be described 
by the theory Ref|3- The reason for positive magnetore- 
sistance is a spin-orbit interaction. In semiconductor het- 
erostructures it is described by the following Hamiltonian 



H{k) = h(T- n{k), 



(1) 



where k is the electron wave vector, er is the vector of 
Pauli matrices, and $7 is an odd function of fc. The spin 
splitting due to spin-orbit interaction Eq. JQ) equals to 
2hn{k). 

In the presence of spin splitting, weak magnetic field 
decreases conductivity. Therefore the effect in such sys- 
tems is called weak antilocalization. Theory of magne- 
toresistance in systems with spin-orbit interaction Eq.lpQ) 
was developed in Ref|3- However the obtained expres- 
sions are valid only for i)weak spin-orbit interaction and 
ii)very low magnetic fields. The first assumption means 
that r2r ^ 1, where r is the scattering time. The second 
condition reads as Zb 3> I, where Ib = ^/h/eB is the 
magnetic length, and I is the mean free path. This so- 
called "diffusion" regime takes place in fields B ^ Btr, 
where 



Bf. 



h 

2^2 



is the "transport" field. 

In high-mobility structures both these conditions fail. 
Due to long scattering times the product fir can be even 
larger than unityi^i^iSii Besides, the transport field is of- 
ten less than 1 mTfS&Sr that is too small range of mag- 
netic fields. This means that particle motion is rather 
ballistic than diffusive. Therefore fitting experimental 
data by the theories Refs0|3 is not always successful.^ 



An attempt to derive the field dependence of anoma- 
lous magnetoresistance for high-mobility structures has 
been performed in Refi6|. However the developed theory 
is correct only for high fields B/Btr ^ {^t)'^ and ignores 
some contributions to the conductivity. 

The aim of the present work is to develop the weak- 
antilocalization theory for systems with strong spin-orbit 
interaction valid for both ballistic and diffusion regimes. 
The magnetic field dependence of the conductivity is cal- 
culated for arbitrary values of B/Btr and Qt, opening a 
possibility to describe anomalous magnetoresistance ex- 
periments and to extract spin-splitting and kinetic pa- 
rameters of high-mobility 2D systems. 



II. THEORY 

There are two fc-linear contributions to the spin-orbit 
interaction Hamiltonian Eq. Q in 2D semiconductor sys- 
tems: the Rashba term and the Dresselhaus term 
rj/j. In heterostructures grown along the direction z \\ 
[001] both vectors flB.D lie in the 2D plane and have the 
following form 



J7/ 



(sinx, -cosx), ^ (cosx, -sinx). (2) 



Here the axes are chosen as x \\ [100], y || [010], and 
tanx = ky/kx. The anomalous magnetoresistance is the 
same if one takes into account the Rashba or the Dres- 
selhaus contribution. Therefore we consider below only 
one term in $7 with an isotropic spin splitting 2TiQ ^ k. 

Retarded and advanced Green functions of a system 
with the spin-orbit interaction Eqs. (0) are 2x2 
matrices in the spin space. In the Landau gauge under 
scattering from a short-range potential, they are given 



G«^^(r,r') = 



Nqs 



ENs±'iTT'/'2.T±ih/2T^' 



(3) 



Here is the Fermi energy, is a number of the Lan- 
dau level, q is the wave vector in the 2D plane, is a 
phase relaxation time, s = ± enumerates two spin states, 
E]Ss is the electron energy, and two-component spinors 
Nqs are the electron wave functions in the presence of 



2 



magnetic field and the spin-orbit interaction Eq. . For 
Rashba spin-splitting, '^Nqs is a superposition of the elec- 
tron states |Af,(7,t) and \N + 1, q, Df^^i.e. with the same 
N + Sz, where Sz is the spin projection onto the growth 
axis. The Dresselhaus spin-orbit interaction mixes the 
states with equal N — Sz- 
In low magnetic fields 



CJc < f^, T ^ < Ep/h, 



(4) 



where uJc is the cyclotron frequency, one can show that 
both magnetic field and spin-orbit interaction result in 
an appearance of phases in the Green functions 

G^-^{r, r') ^ G^'^iR) exp [i(^(r, r') - Wa ■ uj{R)]. (5) 

Here R= {X,Y) =r- r', G^'^(i?) are the Green func- 
tions at i? = and = 0, and (p(r, r') = (x + x'){y' — 
y)/2l'^. The vectors are determined by the sym- 

metry of a spin-orbit interactionii 

u;n{R) = ^{Y,-X), uj^{R) = ^{X-Y). (6) 

The weak-localization correction to conductivity is de- 
termined by interference of paths passing by a scatter- 
ing particle in opposite directions. The amplitude of 
this interference, Cooperon, depends on four spin indices: 
Ca/3,7(5(^, f"')- Here a and /3 (7 and 5) are the spin states 
of a particle before and after passing the path between 
the points r and r' in the 2D plane forward (backward). 
The Cooperon satisfies the matrix equation 

C(r,r') = — P(r,r')+ / driP(r, ri)C(ri, r'), (7) 
rriT J 

where m is the electron effective mass and 



-Gf,(r,r')G4(r,r' 



The expansion coefBcients are given by 



P{N, N') 



j dr Po{r) exp [-2irS • w(r)]Fjvjv 



where 



■N 



Here t = r/ls, tan0 = y/x, and LfJ are the associated 
Laguerre polynomials. At = 0, P{N,N') ~ Snn'- 
Finite spin splitting leads to nonzero values of P{N, N') 
with \N-N'\ < 2. 

Expanding the Cooperon ^C{r,r') in series © as 
well, we obtain the following infinite system of linear 
equations for its expansion coefficients 

C{N,N') = P{N,N') + ^P{N,Ni)C{Ni,N'). (10) 

In order to solve this system, we turn to the represen- 
tation of total angular momentum of two particles S: 
a7 Siris, where S = 0, 1 is the absolute value of S, 
and rris is its projection onto the z axis (|ms| < S). The 
pair of particles with S* = is in the singlet state while 
5 = 1 corresponds to the triplet one. 

Spin-orbit interaction Eq. with only one contri- 
bution Eq. remains the energy spectrum isotropic. 
Therefore the particles being in the singlet and triplet 
states do not interfere, and there are two uncoupled 
Cooperons corresponding to the triplet and singlet, Ct 
and Cs- 

The singlet part is independent of the spin splitting 
and can be found from Eq. H10() as for = 0: 



Cs{N,N') 



P, 



N 



1-Pn 



JNN' 



(11) 



where 



is the probability for an electron to propagate from r to 
r' forward and backward.^ It follows from Eq. Q that 

P{r,r') ^ Pq{R) exp [2i(/7(r, r') -21x5 •a;(fi)], (8) 

where Sa^^ps = {(Tap + (Tjs)/2 is an operator of the total 
angular momentum of two interfering particles, and 



Po{R) 



exp i-R/l) 
2ttRI 



is the value of P in the absence of a magnetic field and 
a spin-orbit interaction. Here the effective scattering 
length I ^ 1/{1 + t/t^). 

In order to find the Cooperon, we expand the matrix 
P into the series over wave functions of a spinless particle 
with the charge 2e in a magnetic field 

P(r,r') = ^ P(iV,iV') <^>r,,{r)<^>%,g{r'). (9) 

NN'q 



P 



N 



I 



dx exp 



Ib 

-X — 

I 



X 

~2 



Ljv(x2). 



The triplet part Ct{N,N') satisfies Eq. ^ with an 
infinite matrix Pt{N^ N'). Ct and Pt are matrices with 
respect to both TV, N' and mg, m'^ = 1, 0, —1. It is crucial 
that Pt can be decomposed into 3x3 blocks. For Rashba 
spin-orbit interaction, it takes place in the basis of the 
states \N,ms) with equal N + ms: |iV- 2, 1), |7V- 1,0), 
|iV, — 1), while for Dresselhaus term this takes place for 
the states with the same N — ms- In both cases the blocks 
in Pt can be obtained by a unitary transformation from 
the following matrix 



P, 



A 



N 



N-2 

R 



-S' 
(1) 

N-2 

-.(2) 



(0) 

Ar-2 



R 



(1) 

N-2 



Pn-1 - 2S 



R 



(1) 

N-l 



(0) 



'(2) 

N-2 

.(1) 

■N-l 



Pn-S 



(0) 
N 



(12) 
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Here 



(„) ^ Ib / Nl 
^ I V (N + my. 



X /" da; exp - y j a;"'L^(x2) sin ( 20ry x 



The triplet part of the Cooperon is expressed via the 
matrix (|12|l as follows: it consists of the blocks Ct{N) 
given by 



(13) 



where / is a 3 x 3 unit matrix. 

The conductivity correction due to weak antilocaliza- 
tion is given by a sum of two terms^ 

a{B) = (Ja+ (Tb, 

where Ca and erf, can be interpreted as backscatter- 
ing and nonbackscattering interference corrections to 
conductivity.^^ They are given by 

= 1^ / / E Ca7,/35(r-, r') (14) 
xJsa{r',r) ■ Jp^{r',r), 



<Jb 



27rmT 



j dr j dr' J dr" ^ C„^,^5(r,r') (15) 



JLir\r'')JUr'',r)Gl{r',r'')GUr'',r) 



+ J|,(r', r") J^^Jr", r)G^;(r', r")G^V^"' ^) 



Appearance of the modified Cooperon C = C ~ fol- 
lows from that only three and more scattering events con- 
tribute to the magnetoconductivity. 
The current vertex is defined as 

J(r-,r')-e j driG^(r, ri)t;(ri)G^(ri, r'), 

where v is the velocity operator in a magnetic field. Sub- 
stituting here the Green functions in the form Eq. 
one can show for low magnetic fields Eq. Q that 

lel 



where = ± iJ'^ , and 6 is an angle between r and 
r'. 

Omitting the rapidly oscillating products G^G^ and 
G'^G"^ and expanding the terms 

K(r,r') = icos9 P{r,r') 

in series Eq. @, we get from Eqs. (|14|I - H16|) the final 
expressions for the conductivity corrections 

(17) 



I 



Air'^h \Ib 



2 OO 



r IZ^Tr ifjvi^wAjv(/-A^v)-Ml8) 



JV=0 



-Tr ^KNKNAN+iil - An+iY 



p. 



N 



P. 



N + l 



l~P, 



N 



l~Pl 



N + l 



The terms with matrices here are the triplet contribu- 
tions which is seen to be of opposite sign in comparison 
to the singlet ones. The matrices K]\[ and K]\[ appearing 
in the expansion of the function K{r, r') are given by 



Qn-2 — >5'^^2 



R 



Kn = 



AO) 

■N-l 
(1) 



(2) 
JV-2 



25* 



(1) 



j(0) 



c.(3) 

(2) 

N-l 



R 



-R 



N 



S 



(1) 

N 



(19) 



K 



N 



where 



^ Qn-2 — S 
-R 

V 



N-2 - '3n-2 
(2) 
JV-2 
c(3) 



-R 



(0) 



Qn-1 — 25* 



'N-2 



-R 



(2) 
N-1 



(1) 

N-l 



;.(!) 

(0) 

■N 



R 



Qn — S 



(1) 

N 



(20) 



Qn — 



1 h 



OO 

X y ^''^p 





Note that the values with negative indices appearing in 
Eqs. lO, (UU, and ^ at iV 0, 1 should be replaced 
by zeros. 

Eqs. p7|l and p8|) yield the weak-antilocalization 
correction to the conductivity in the whole range of 
classically-weak magnetic fields and for arbitrary values 

of nr. 



III. LIMITING CASES 



In the hmit of zero spin splitting, 5^'^ — — 0, 



"^a/3('"' ^ ^ ~ p. ^ [^0/3(^1 r ) + G^p{r, r )] , (16) ^}jg matrices Ajy, ifjv, and Kn became diagonal, and we 
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obtain 



^_ (]_ 



2 oo 



p3 



^ 1 



R / ' • 1 — _P/V 



(21) 



Oh = 



2 oo 



7V=0 



TV 



1 -P/ 





Eqs. I|21|l and l|22|l coincide with the results of non- 
diffusive theory developed for = in Ref^gj. 

In the diffusion regime, when B <§C one can calcu- 
late the difference between the conductivity in the pres- 
ence and in the absence of magnetic field, ^a{B). Mak- 
ing use of standard approximations valid in the diffusion 
regime, we obtain from Eq. H17() 



(22) The matrices here are given by 



(0) 



R 



(1) 



(1) 

X 

(2) 



(0) 



(2) 

X 

(1) 



(28) 



R 



(1) 



Px Si 



(0) 



AfTd.//(B) 









(i 


^x Sx 


R^' ^ 




FsiB) 


(23) 


Kx = 






Qx '^Sx 










[ 


c(l) 
^x 




Qx "S": 



(29) 



Here the singlet contribution is given by 

Fs{B)^^>{l/2 + b^)-\nb^, 



(24) 



where ^' is the digamma-function. The expression for the 
triplet term is as foUowsi^ 



Ft{B) 



J_ _ 2ao + 1 + bs 
ao ai{an + bs) -2C,bs 




(e + 2)a^ + 2aNbs - ^ - 2(2jV + l%bs 
{aN + bs)aN-iaN+i - 2Cbs[{2N + l)ajv - 1] 



e + 3 



In I 



e+1 



2 ^ " 2 

Here iV + 1/2 + 6^ + fe^, 

T 2(r2T)2 



In (6^ -f 26,). (25) 



^ B ' ' 1 + (2f7r)2 B ' 

e = [l + 2(^^T)2]-^ C= [l + (2^^T)2]-^ 

Equations (|23|I - H25|) generalize the diffusion theory Ref0 
to the case of arbitrary strong spin-orbit interaction. In 
the limit Qt <ti 1, these expressions pass into the results 
of Refj2r 

The conductivity correction in zero magnetic field, 
(t(0), can be obtained from Eqs. (|17|l and p8|) by pass- 
ing from summation over N to integration and using the 
following asymptotic valid for ^ 1 

1 



=a;'"L^(a;2) « Jmi2xVN). 



As a result, we get for ^ t 



(26) 



|i JdxTr[A^(/-A.)-i]-ln^|, 
' 



Kx 
where 

Px = 





^x Sx 


p(0) 

— nx 








— ^x 


Qx 2Sx 




(30) 


V 


c(3) 

^X 


p(2) 

— nx 


Qx Sx j 





^(l + r/r^)2+x 

OO 

y dyexp(-j/)J„(yV^)sin2(r!r2/), 



1 

V2 



exp {-y)Jrn{yVx) sin (2riTj/). 



In Refs,15 cr(0) has been analyzed in the diffusion ap- 
proximation In (r0 /r) 3> 1 which is hardly realized prac- 
tically. 

In a magnetic field B ^ {flT)^Btr, the conductivity be- 
comes independent of fl. The reason is that in so strong 
field the dephasing length due to magnetic field ~ Ib 
is smaller than one due to spin-orbit interaction, I /fir. 
As a result, the particle spins keep safe at characteristic 
trajectories. The conductivity for any finite Qt has the 
zero-r2 asymptotic Eqs. (|2l}, (1221) • For fir < 1 this de- 
pendence is achieved at _B ^ Btr- In high magnetic field 
B ^ Btr, {^T)'^Btr, the conductivity correction has the 
high-field asymptotio^ 



'^hf{B) 
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FIG. 1: Conductivity correction (solid curve) at fir — 1, 
t/t^ = 0.01. Dashed curves represent the backscattering (tTa) 
and nonbackscattering (at) contributions, dotted curves show 
the results of diffusion and high-field approximations. 



IV. RESULTS AND DISCUSSION 

In Fig. n] weak-antilocalization correction to the con- 
ductivity for rtr = 1 is shown by a solid line. In low 
fields the conductivity decreases and reaches a minimum 
at some B = Bmin- Then the field dependence asymp- 
totically tends to zero. Dashed curves in Fig. ^ represent 
the backscattering (ctq) and nonbackscattering (df,) con- 
tributions. One can see that df, can reach almost 25% of 
\<7a\, therefore the nonbackscattering correction should 
be taken into account when fitting experimental data. 
The dotted curves in Fig. ^ show results of the diffusion 
and high- field approximations.^^ One can see that the 
former is valid in a narrow region B < 0.5 Btr- The 
high-field asymptotic holds true only for B > 100 Btr- 
This proves importance of non-diffusion theory for high- 
mobility structures. 

In Fig. 121 the conductivity correction is plotted for dif- 
ferent strengths of spin-orbit interaction. One can see 
that for fir ^ 1, in accordance with results of the previ- 
ous Section, a{B) coincides with the zero-il dependence 
for B > Bmin- The asymptotic ahf{B) is reached at 
B « 100 Btr for all finite values of fir. The positions of 
minima in the curves are shown in the inset. One can see 
that Bmin almost linearly depends on the spin splitting 
at fir > 0.8. Fitting yields the following approximate 
law 



Br, 



(3.9 VlT - 2)Btr 



In the limit r2r — > oo, the triplet state with rus = 
does not contribute to the conductivity. The correspond- 



FIG. 2: Conductivity correction for different strengths of spin- 
orbit interaction at t/t^ = 0.01. The inset represents the 
positions of minima in the magnetoconductivity. 



ing dependence is presented in Fig. |21 One can see a de- 
crease of conductivity in the whole range of magnetic 
fields. At B ^ Btr, the correction tends to zero as 
0.035 eV?i y/Btr/B. 

In experiments, the difference Aa{B) = a{B) — a(0) 
is measured. Since (7(B) tends to zero at B oo for 
any fir, one can extract cr(0) from the saturation value 
of Act at _B ^ oo. In Fig. O the zero-field value of the 
conductivity correction a{0) is plotted as a function of 

riT. 

Fig. 121 shows how spin-orbit interaction changes the 
sign of weak-localization correction to conductivity. At 
fir — 0, all three triplet states and a singlet one yield 
contributions of the same absolute value. As a result, 
the zero-field correction is given by 



■In^, 



■In 2. 



In the opposite limit Qt — > oo, the triplet contribution is 
partially suppressed by the spin-orbit interaction. Cal- 
culation shows that the corrections reach the following 
values 



= (0) 



(0.57 + In ^) 



^(0) 



-0.43- 



(31) 

One can see that ct(0) changes its sign and reduces its 
absolute value when Qt increases from zero to infinity. 

It follows from Fig. that the magnetoconductivity 
Af7(i?) is an alternating function at small Or, while at 
large values of the spin splitting Aa{B) is negative in the 
whole range of classically weak magnetic fields. 



6 



as 

u 

u 
o 



I 

o 



0.1 



0.0 



-0.1 



-0.2 



1 1 


1 






// 


- - . 


// 


- 


// 




^/ (7 

^ a 







0.01 



0.1 



10 



100 



FIG. 3: Zero-field correction to the conductivity for r/r^ = 
0.01 (solid curve). Dashed curves represent the backscattering 
(oa) and nonbackscattering [ph) contributions. 



V. CONCLUSION 

In the present paper the anomalous magnetoresistance 
is calculated for 2D systems with only Rashba or only 
Dresselhaus spin-orbit interaction. In both cases the 
spin splitting is isotropic in fc-space and characterized 
by one constant In the presence of both types of spin- 
orbit interaction, Eqs.®, and lf7|l- lfTT|) hold true with 
u; = LOfi -|- ujy). However Pt{N,N') is not divided into 
finite blocks as Ea. H12() . and one should use an infinite 
matrix for calculation of the triplet contribution to the 
conductivity in this case. 

The problem has an analytic solution if the Rashba 
and Dresselhaus spin splittings are equal to each other. 
In this case the magnetoconductivity is positive in the 
whole range of magnetic fields like in systems without 
spin-orbit interaction. This result has been previously 
obtained in the diffusion approximation for B <C Btr^ 



For magnetic field of arbitrary strength, the dependence 
(7{B) is given by Eas. (|^ . if^ . This can be proved by 
noting that at H-r = flo the vector J7(fc) is directed along 
the same axis for all fc. As a result, the energy spectrum 
consists of two identical paraboloids shifted relative to 
each other in the direction of Both these spin sub- 
bands independently yield equal conductivity corrections 
coinciding with those for spinless case. The same result 
takes place for symmetrical [110]- and [113]-grown quan- 
tum wells. 

Application of an in-plane magnetic field B|| destroys 
weak antilocalization. It has been demonstrated experi- 
mentally that the magnetoconductivity minimum disap- 
pears in the presence of B^\M^ Weak antilocalization 
in a tilted magnetic field can be also described by the 
present theory. Parallel field influences the anomalous 
magnetoresistance due to two microscopic reasons. First, 
an in-plane field results in additional dephasing due to or- 
bital effects. ^^-^^ They can be taken into account as 
dependent corrections to t^. Second, an in-plane field 
induces finite Zeeman splitting. This results in a mixing 
of the singlet and triplet states, which makes the ma- 
trix P{N, N') in Ea. (|10ll infinite. However if the Zeeman 
splitting is much smaller than hfl, then By affects only 
the singlet state. It leads to another correction to 
which should be taken into account only in Cs, Ea. (|ll|l . 
Both the Zeeman and the orbital corrections to the de- 
phasing rate can be extracted from the fit of experimental 
data by Eas. (|17|l . 118(1 . Inclusion of the dephasing correc- 
tions into Eas. (|26|l . (|27|l allows one to describe anomalous 
magnetoresistance in pure in-plane field as well. 

In conclusion, the theory of weak antilocalization is de- 
veloped for high-mobility 2D systems. Anomalous mag- 
netoconductivity is calculated in the whole range of clas- 
sically weak fields and for arbitrary values of spin-orbit 
splitting. 
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